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Experimental and numerical data for the heat transfer as a function of the Rayleigh-, Prandtl-, 
and Rossby numbers in turbulent rotating Rayleigh-Benard convection are presented. For relatively 
small Ra ~ 10^ and large Pr modest rotation can enhance the heat transfer by up to 30%. At 
larger Ra there is less heat-transfer enhancement, and at small Pr < 0.7 there is no heat-transfer 
enhancement at all. We suggest that the small-Pr behavior is due to the breakdown of the heat- 
transfer-enhancing Ekman pumping because of larger thermal diffusion. 
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Turbulent convection of a fluid contained between 
two parallel plates and heated from below, known as 
Rayleigh-Benard convection (RBC), continues to be a 
topic of intense research P, 0, Q . A particularly interest- 
ing variation of RBC is the case where the sample is ro- 
tated about a vertical axis at an angular speed Vt. That 
system is relevant to numerous astro- and geo-physical 
phenomena, including convection in the arctic ocean 
in Earth's outer core [5], in the interior of gaseous giant 
planets [6|, and in the outer layer of the Sun [7|. Thus 
the problem is of interest in a wide range of sciences, 
including geology, oceanography, climatology, and astro- 
physics. 

It is widely understood [8] that rotation tends to sup- 
press convective flow, and with it convective heat trans- 
port, when the rate of rotation is sufficiently large. How- 
ever, at modest rotation rates, experiments [J^ [13, . 
and numerical simulations [l2|, [l^, [lH have shown 
that under certain conditions the heat transport can also 
be enhanced^ before it rapidly decreases for stronger ro- 
tation. This enhancement has been ascribed to Ekman 
pumping 1^, l^^IS]* Due to the rotation, rising or 
falling plumes of hot or cold fluid are stretched into ver- 
tical vortices that suck fluid out of the thermal boundary 
layers (BL) adjacent to the bottom and top plates. This 
process contributes to the vertical heat flux. There are, 
however, some experiments [l9j and numerical simula- 
tions [ij] that show hardly any enhancement of the heat 
flux at modest rotation rates. 

In the present paper we determine systematically as a 
function of the Rayleigh number i?a, Prandtl number Pr, 
and Rossby number Ro (all to be defined below) where 
the heat-flux enhancement occurs. We present both ex- 
perimental measurements and results from direct numer- 



ical simulation (DNS). They cover different but overlap- 
ping parameter ranges and thus complement each other. 
Where they overlap they agree very well. We find that 
in certain regimes the heat-flux enhancement can be as 
large as 30%; this raises the possibility of relevance in 
industrial processes. Even more remarkably, we observe 
a heretofore unanticipated strong dependence of this en- 
hancement on Pr as well as on Ra. 

For given aspect ratio F = D/L {D is the cell di- 
ameter and L its height) and given geometry (and here 
we will only consider cylindrical samples with F = 1), 
the nature of RBC is determined by the Rayleigh num- 
ber Ra = l3gAL^ /{f<iu) and by the Prandtl number 
Pr = u/tv. Here, P is the thermal expansion coefficient, g 
the gravitational acceleration, A = T5 — the difference 
between the imposed temperatures T5 and Tt at the bot- 
tom and the top of the sample, respectively, and u and 

the kinematic viscosity and the thermal diffusivity, re- 
spectively. The rotation rate ^} (given in rad/s) is used 
in the form of the Rossby number Ro = (3gA/L/{2Q). 

The convection apparatus was described in detail as 
the "medium sample" in Ref. |20|. Since the previous 
measurements [21] it had been completely dis- and re- 
assembled. It had copper top and bottom plates, and 
a new plexiglas side wall of thickness 0.32 cm was in- 
stalled for the current project. The sample had a diam- 
eter D = 24.8 cm and a height L = 24.8 cm, yielding 
F = 1.00. The apparatus was mounted on a rotating 
table. We used rotation rates up to 0.3 Hz. Thus the 
Froude number Fr = Q?{L/2)/g did not exceed 0.05, 
indicating that centrifugal effects were small. Cooling 
water for the top plate and electrical leads were brought 
into and out of the rotating frame using Dynamic Sealing 
Technologies feed-throughs mounted on the rotation axis 
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FIG. 1: The ratio of the Nusselt number Nu(Q) in the pres- 
ence of rotation to Nu{Q = 0) for Pr = 4.38 {Tm = 40.00° C). 
(a): Results as a function of the rotation rate in rad/sec. (b): 
The same results as a function of the Rossby number Ro on a 
logarithmic scale. Red solid circles: Ra — 5.6x10^ (A = 1.00 
K). Black open circles: Ra = 1.2 x 10^ (A = 2.00 K). Pur- 
ple solid squares: Ra = 2.2 x 10^ (A = 4.00 K). Blue open 
squares: i^a = 8.9x10^ (A = 16.00 K). Green solid diamonds: 
Ra = 1.8 X 10^° (A = 32.00 K) 



above the apparatus. All measurements were made at 
constant imposed A and 1], and fluid properties were 
evaluated at = {Tt + T}))/2. Data were taken at 
60.00°C (Pr = 3.05, ty = L'^/u = 1.27 x 10^ sec), 
40.00°C (Pr = 4.38, ty = 9.19 x 10^ sec), 24.00°C 
(Pr = 6.26, ty = 6.69 x 10^ sec) and 23.00°C (Pr = 
6.41, ty = 6.56 X 10^ sec). In a typical run the system 
was allowed to equilibrate for three or four hours, and 
temperatures and heat currents were then averaged over 
an additional three or four hours and used to calculate 
Ra and the Nusselt number Nu = QL/{XA) {Q is the 
heat-current density). 

Measurements of the Nusselt number without rota- 
tion, Nu{n = 0), over the range 5 x 10^ < < 10^° 
(1 ^ A ^ 20K) agreed within estimated systematic er- 
rors of about 1% with previous results |21|] obtained in 
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FIG. 2: The ratio Nu{^)lNu(Sl = 0) for Ra = 1.2 x 10^ 
as function of Ro on a logarithmic scale. Red solid circles: 
Pr = 3.05 {Tm = 60.00° C). Black open circles: Pr = 4.38 
(Tm = 40.00° C). Purple solid squares: Pr = 6.41 {Tm = 
23.00° C). Blue stars: Numerical simulation for Ra = 1.0 x 10^ 
and Pr = 6.4 from Ref. 



the same apparatus. The ratio Nu{ft)/Nu{ft = 0) is 
shown in Fig. [1^ as a function of the rotation rate Q. 
Those results are for = 40.00°C, where Pr = 4.38. 
The enhancement of Nu due to modest rotation is clearly 
seen at all Ra. It is larger at the smaller Ra. 

In Fig. [1)3 we show the same data as a function of the 
Rossby number Ro. At large Ro (small the data must 
approach unity, and indeed they do. As Ro decreases, 
Nu is first enhanced, but then reaches a maximum and 
decreases as expected. The maximum of Nu{Q) /Nu{Q = 
0) occurs at larger Ro for larger Ra, and the value of this 
ratio at the maximum diminishes with increasing Ra. 

In Fig. [2] the results at constant Ra are shown as a 
function of Ro for several values of Pr. Also shown are 
the DNS results from Ref. [15]; these data agree rather 
well with the experimental data at nearly the same Pr 
and Ra. One sees that the enhancement of Nu at large 
Ro is nearly independent of Pr; but as Ro decreases be- 
low unity a strong Pr dependence manifests itself. As Ro 
decreases, the depression of Nu sets in earlier for smaller 
Pr, and the maximal relative heat transfer enhancement 
is smaller. 

In the DNS we solved the three-dimensional Navier- 
Stokes equations within the Boussinesq approximation. 
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with V • u = 0. Here £ is the unit vector pointing in 
the opposite direction to gravity, u the velocity vector, 
and the non-dimensional temperature, < ^ < 1. 



3 




FIG. 3: The ratio Nu{Q)/Nu{Q = 0) as function of Ro on 
a logarithmic scale. Red solid circles: Ra — 2.73 x 10^ and 
Pr = 6.26 (experiment). Black open circles: Ra = 2.73 x 10^ 
and Pr = 6.26 (DNS). Blue solid squares: Ra = 1 x 10^ and 
Pr = 6.4 (DNS) O. Red open squares: Ra = 1 x 10^ and 
Pr = 6.4 (DNS). Green open diamonds: Ra ^ 1 x 10^ and 
Pr = 20 (DNS). Black stars: Ra = 1 x 10^ and Pr = 0.7 
(DNS). 



Finally, P is the reduced pressure (separated from its 
hydrostatic contribution, but containing the centripetal 



contributions): P = p — r / {8Ro with r the distance 
to the rotation axis. The equations have been made non- 
dimensional by using, next to L and A, the free-fall ve- 
locity U = \^l3gAL. The simulations were performed 
on a grid of 129 x 257 x 257 nodes, respectively, in the 
radial, azimuthal and vertical directions, allowing for a 
sufficient resolution of the small scales both inside the 
bulk of turbulence and in the BLs (where the grid-point 
density has been enhanced) for the parameters employed 
here 14, T^. Nu is calculated as in Refs. and its 

statistical convergence has been controlled. 

The numerical results for Nu as function of Ro for sev- 
eral Ra and Pr are shown in Fig. [3l For Ra = 2.73 x 10^ 
and Pr = 6.26 experimental data (not previously shown 
in Figs. [Hand [2]) are plotted as well, showing near-perfect 
agreement with the numerical results. This gives us con- 
fidence that the partial neglect of centrifugal effects in 
the simulations [namely, neglecting the density depen- 
dence of the centripetal forces, which in the Boussinesq 
equations show up as — 2Fr rOr [22[ (with the radial unit 



vector r )] is justified, as already found in [l3|, ll5|, be- 
cause Fr <C 1 . Neither do the experimentally unavoidable 
finite conductivity of the top and bottom plates 13, [1^ 
and the side- wall conductivity 0, seem to matter in 



Ref. [3^ for the non-rotating case. 

As was seen for the experimental data in Figs. [1] and 
[2 the numerical results in Fig. [3] also reveal a drastic 
dependence of the Nusselt-number enhancement on Pr. 
For large Pr > 6 the enhancement can be as large as 
30% for Ra ^ 10^ and Ro ^ 0.1. However, there IS no 
enhancement at all for small Pr < 0.7. This trend is 
further elucidated in Fig. [H where we show the Nusselt- 
number enhancement as function of Pr for three different 
Ro and Ra = 10^. 




this regime of parameter space, as already explained in 



FIG. 4: Numerical result for the ratio Nu{Q)/Nu{0) as func- 
tion of Pr for Ra = 10^ and Ro =1.0 (red open diamonds), 
Ro = 0.3 (black open circles), and Ro = 0.1 (blue open 
squares) . 



Why does the heat-transfer enhancement through Ek- 
man pumping at modest rotation rates break down as Ro 
decreases below a Pr-dependent typical value, as seen 
in Figs. [1] to [3l? To obtain a hint, we visualized (see 
Fig. [5]) the three-dimensional temperature iso-surfaces 
for Pr = 0.7 and for Pr = 6.4 at both Ro = 0.30 and 
Ro = 00, at Ra = 10^. While for the larger Pr = 6.4 
case the temperature iso-surfaces organize to reveal long 
vertical vortices as suggested by the Ekman-pumping pic- 
ture, these structures are much shorter and broadened for 
the low Pr = 0.7 case, due to the larger thermal diffu- 
sion which makes the Ekman pumping inefficient. This 
would imply enhanced horizontal heat transport which 
should also lead to a steeper gradient of the mean tem- 
perature in the bulk. Indeed, in the DNS, we find that 
when Ro becomes small enough, the bulk of the fluid 
displays an increasingly destabilizing mean temperature 
gradient (see Fig. [6|), which of course must be accompa- 
nied by a reduction of the mean temperature drop over 
the thermal BLs and thus a Nusselt-number reduction. 
The first manifestation of the enhancement of the mean 
destabilizing vertical temperature gradient agrees with 
the onset of relative heat-transfer reduction in figure [31 
and thus supports this explanation. 

Along the same line of arguments one may also expect 
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FIG. 5: 3D visualization of the temperature isosurfaces at 
0.65 A (red) and 0.35 A (blue), respectively, for Pr = 0.7 
(upper plot) and Pr = 6.4 (lower plot) for Ra = 10^ and 
Ro = oc(left) Ro = 0.30 (right). 
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FIG. 6: The horizontally averaged vertical temperature gra- 
dient dO/dz at the sample mid-plane as a function of Ro for 
Ra = 1 X 10^. Blue open squares: Pr — 0.7. Red open circles: 
Pr — 6.4. Black open diamonds: Pr — 20. 

that the Ra dependence of the reduction of Nu at small 
Ro seen in Fig. [1] is attributable to relatively less efficient 
Ekman pumping at higher Ra\ The enhanced turbulence 
may lead to a larger eddy thermal diffusivity, promoting 
a homogeneous mean temperature in the bulk. Again 
this would make Ekman pumping relatively less efficient 
and reduce the peak in the relative Nusselt number. 

Another interesting aspect of our data is that within 
our experimental or numerical resolution there is no heat- 



fiux enhancement for Ro > 2 for any Ra or Pr. As 
already noticeable from the data of Ref. [lH, the heat- 
flux enhancement first becomes resolved as Ro decreases 
below about two. Future work has to reveal whether 
this onset coincides with some reorganization of the flow. 
Similarly, it remains to be analyzed how the dramatic de- 
pendence of Nu on Ra and Pr for modest Ro is reflected 
in the flow fields. 
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